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"j^ . Abstract 

Kinetic aspect of Landauer's blowtorch effect is investigated for a model double- 
well potential with localized heating. Using the super symmetric approach, we de- 
rive an approximate analytical expression for the equilibration rate as function 
of the strength, width and the position of the hot zone, and the barrier height. 
We find that the presence of the hot zone enhances the equilibration rate, which 
is found to be an increasing function of the strength and width of the hot zone. 
Our calculations also reveal an intriguing result, namely, that placing the hot 
zone away from the top of the potential barrier enhances the rate more than when 
it is placed close to it. A physically plausible explanation for this is attempted. 
The above analytical results are borne out by detailed numerical solution of the 
associated Smoluchowski equation for the inhomogeneous medium. 



1 Introduction 

In a now influential paper |l[ on the relative stability, i.e., relative occupation, of the com- 
peting local energy minima for a system far from equilibrium, Landauer pointed out the 
globally determining role played by the non-equilibrium kinetics of the unstable intermedi- 
ate states even as these are very rarely populated. More specifically, for the case of two 
local energy minima, i.e., a bistable potential, he showed that the application of a localized 
heating at a point on the reaction coordinate lying between the lower energy minimum and 
the potential barrier maximum can raise the relative population of the higher lying energy 
minimum over that given by the thermal Boltzmann factor exp(—A E/ksT). This is the 
so called 'blowtorch' effect |J associated with a nonuniform thermal bath. It generalizes 



the problem of escape of a Brownian particle over a potential barrier under the influence of 
equilibrium thermal fluctuations, studied originally by Kramers [2-4], to the case of nonuni- 
form temperature along the reaction coordinate. Later, in a related context, it has been 
shown that a state-dependent diffusion coefficient can produce maxima in the probability 
distribution at points which are not the potential minima [5-7]. 



A decade later, van Kampen || derived an equation appropriate for the description of dif- 
fusion in an inhomogeneous medium where he developed a stochastic treatment for the case 
of non-uniform temperature. Apart from justifying Landauer's conjecture, he also showed 
that there could be a net current when the particles are allowed to diffuse back through 
an alternate route bypassing the hot zone. The problem has also been treated by Biittiker 
around the same time ||. He showed that a net current is possible, even in the absence 
of an externally applied field, provided both the potential and the state-dependent diffu- 
sion constant are periodic with a relative phase difference. Later Landauer [|TUJ considered 
this aspect of the problem again in the light of van Kampen's work. Sinha and Moss 
have verified Landauer's conjecture by computer simulation. This has also been applied to 
thermal activation in a superconducting ring with a weak link where transitions produce 
temperature changes |i~2[ . Indeed, a whole new field of research has emerged centering on 
the idea of the possibility of directed motion out of noisy states under diverse athermal 
driving conditions, generically subsumed under 'thermal ratchets' [13-18], and traceable to 
the original 'blowtorch theorem' of Landauer in the sense that the latter may be viewed as 



injection of noise [19 



Most of the investigations based on the blowtorch effect study the influence of space- 
dependent temperature on the steady-state relative occupations of the energy minima. To 
the best of authors' knowledge, there has been no analytical attempt to study the kinetic as- 
pect of the system, specifically, the calculation of the longest of relaxation time in a bistable 
potential in the presence blowtorch (hot zone). In this work, we address this problem and 
calculate the equilibration rate for a simple bistable potential with the localized heating posi- 
tioned somewhere in the (unstable) intermediate region. The associated Smoluchowski equa- 
tion describing the physical situation will be dealt with using the supersymmetric (SUSY) 
approach [20-22]. In doing so, we have adapted the supersymmetric potential approach 
normally applicable to the original Kramers barrier-crossing problem, to the case when the 
temperature is nonuniform along the reaction coordinate. 



The rest of the paper is organized as follows. Section 2 contains a brief introduction to the 
SUSY method for extracting the longest relaxation rate (i.e., equilibration rate) from the 
Smoluchowski equation for the case where temperature is space-dependent. In section 3, 
applying the SUSY method to a model double-well potential with a hot zone, we obtain a 
matrix equation determining the equilibration rate. In section 4, we present an approximate 



analytical expression obtained from the matrix equation. The latter is also solved numerically 
for the sake of comparison. As an independent check for these results, we have solved the 
Smoluchowski equation numerically and obtained the equilibration rate. The results we 
obtained from the above three ways for the equilibration rate as a function of different 
parameters characterizing the hot zone are, then, discussed. We devote the last section to 
summary and conclusions. 



2 The methodology 



The Smoluchowski equation describing the kinetics of a Brownian particle in an inhomoge- 
neous medium is given by |§ 



dP(x,t) 
dt 



d_ 

dx 



U'{x)P{x,t) + ^-(T{x)P{x,t)) 



where the mobility, /i, and temperature, T, are in general space-dependent. (We take Boltz- 
mann's constant, ks = 1.) Here, P(x,t) is the probability density of finding the particle at 
position x at time t and U(x) describes the potential profile and the prime on U denotes the 
derivative with respect to x. 



Some remarks on the use of the Smoluchowski equation, Eq. ([]]), for the case of inhomoge- 
neous medium are in order at this point as the latter has been a subject of some debate 
that still continues ^ [10], [19|, |24|. Landauer has argued J[0J [19] for a generalization of 



the Smoluchowski equation where the second term on the RHS of Eq.([l]) is to be replaced 
by P{x,t) {^f)+aT{x) (^*). The parameter a was shown to depend on the physical 
conditions to be imposed across the temperature discontinuity for no net current. The pa- 
rameter a = 0.5 corresponds to a particle interacting with the thermal bath for which the 
particle velocity is taken to be proportional to \T or P(x) oc -^==. O n the other hand, 

a = 1 corresponds to the case when particles equilibrate via collision (pressure equilibra- 
tion), i.e., P(x) oc yT^y. However, a direct derivation based on phase space Smoluchowski 
equation by van Kampen gave Eq.([l]) corresponding to a — 1. This is also supported by the 
work of Jayannavar and Mahato |25[ based on a microscopic treatment of the thermal bath 
as a set of harmonic oscillators. In our work, we will continue to use Eq.([l]) as providing a 
physically valid description of the problem under consideration. 



We consider the case where temperature is space-dependent but the medium is homogeneous 
so that the mobility, fi, is taken to be constant. Similar assumption of constant fi has been 



taken by Sinha and Moss in their work pT |. Then, the corresponding Smoluchowski equation 
is 



dP(x,t) d 



dt 



^ dx 



U'(x)P(x,t) 



A(T(*)P(M)) 



(2) 



The stationary solution of Eq. (pj) is given by 

C 



Pss(x) 



T(x 



-exp 



— oo Tlx) 



(3) 



where C is the normalization constant. In this work, we choose a simple temperature profile 
whose value is constant both outside and inside the hot zone with T = Tq and T = To + T& 
respectively, where T b is the excess temperature above the constant background value T . 
(See Eq.(|l(]).) Using the following transformation 



POM) 



00) 

T(x) 



-exp 



U'(x) 
2T(x) 



dx 



-xt 



we convert Eq.(H) to a Euclidean Schrodinger equation 



H 4 



A + A(f) + = E+cj) + . 



(4) 



(5) 



Here, E A 



jTjr and the operators A and A + are given by 



A 



fd U'(x) 
dx 2/T 



A+ = _fd U^x) 
dx 2/T ' 



(6) 



In Eq.(||), the parameter / reflects the excess temperature defined by / = \/l + s, with 
Thus, for the interval of x outside the hot zone, the corresponding operators have 



To' 



/ = 1. The Hamiltonian H + corresponds to the motion of a particle in the potential 

(U'(x)\ 2 U"(x) 



VJx) 




2T(x) 



(7) 



For a high barrier, the longest relaxation rate is determined by the smallest nonzero eigen- 
value Ai associated with the first excited state of Eq.((^). On the other hand, this eigenstate 
is degenerate with the ground state 4>°_ of the 'supersymmetric partner potential' V_ (x) given 
by 



VJx) 




+ 



U"(x) 
2T(x) 



such that 



H 



•.o 



AA + <f) 



o 



•.o 



(9) 



where E. 



-^r. Thus, the problem of finding the equilibration rate amounts to finding 
the ground state eigenvalue of this 'partner' potential. It is worth pointing out here that 
for the temperature profile we use, we exclude the two points where T(x) is discontinuous. 
However, following the line of attack used for obtaining the solution of the Schrondinger 
equation in the presence of a delta function potential, we devise a technique of relating the 
wave functions on either sides of these points (Appendix A). 



Here we would like to mention that the SUSY method is a powerful analytical method used 
for simple model potentials such as the W-potentials. The applicability for calculating the 
Kramers' escape rate in a W-potential for a uniform temperature case has been demonstrated 
by Schonhammer [22]. We found this method to be particularly useful while investigating the 
influence of barrier subdivision on the escape rate where we have demonstrated the existence 



of optimal number of subdivisions that maximizes the escape rate [26 



3 The model and its solution 

We consider a simple bistable potential in the form of a symmetric W-potential which is 
piecewise linear having the same magnitude in slope. It is described by the barrier height Uq 
and the distance 2L between the two potential minima located at x — ±L on either side 
of the origin (Fig. la). The localized hot zone of a certain width is taken to be positioned 
somewhere between the two potential minima. We assume a simple temperature profile of 
the following form for the heat bath: 



T{x) =T + T b 



6(x-4 + y)-e(x-4- y) 



(10) 



Here, Q(x) is the Heaviside function, T is the background (constant) temperature and T& is 
the excess temperature of the hot zone. The parameters db and w b specify, respectively, the 
distance of the mid-point of the hot zone from the barrier top and the width of the hot zone 
(Fig. lb). 



The corresponding supersymmetric 'partner' potential VL(a?) with the hot zone is given by 
V.{x) = ^ [6{x + L ) - S(x) + S(x - Lo)] + ( . (11) 

This relation holds good at all points except at the two points of temperature discontinuity. 
Here, u = |jr and u\ = j^. It is clear that the potential V-(x) takes a constant value 

(l^) outside the hot zone and (j^) within. Two additional repulsive and one attractive 
delta function potentials are superposed at different locations as shown in Fig. lc. Changing 
the variable x to y = gives a new dimensionless Hamiltonian h- = LqH_ such that 

M-(V) = e 6 0° (y) (12) 

whose potential is 

V-(y) = 2u [6(y + 1) - S(y) + 5(y - 1)] + (w ,i) 2 (13) 

Here, is a dimensionless quantity equal to L^E_ or equivalently -jjjr-- In the scaled form, 
the position of the mid-point of the hot zone from the top of the barrier and the width of 
the hot zone are defined by d — jf- and w — j^, respectively (see Fig. lb). Henceforth, we 



use only the reduced variables. Transfer matrix method has been used to find the equation 
governing the eigenvalue e&. Below, we outline the procedure leaving the details to the Ap- 
pendix A. 



The ground state wave function, 0° , is taken to have the form A n e~ k( - y ~ yn ^ + B n e k{:!l ~ yn) 
located with respect to the positions y = y n of each of the delta function potentials and at 
points of discontinuity in the temperature profile. The amplitudes A n , B n and A' n , B' n on 
either sides of these special points are shown in Fig. lc. For regions outside the hot zone, 
the wave vector k is given by 

k = \jul- e b , (14) 
while for within the hot zone it is given by 



h = p-rr s < 15) 

The amplitudes A' A and B' A that appear just on the left side of the positive delta function 
potential located at y = —1 are related to A and B that appear on the right side of the 
positive delta function potential located at y — 1 through a transfer matrix M such that 

The matrix M is the product of successive transfer matrices which are derived in the Ap- 
pendix A. Since we are looking for a bound state solution, A' 4 = B = 0. This implies that 
the 11-matrix element of M must be zero, i.e., 

M u = 0, (17) 

the solution of which gives the value of e& and, hence, the equilibration rate. 



4 Results and discussion 

The effect of the hot zone on the equilibration rate for the given potential is studied in terms 
of the three parameters characterizing the hot zone: (i) the relative degree of hotness of 
the hot zone with respect to the rest of the heat bath, s — ^, (ii) the scaled width of the 
hot zone, w = and (iii) the scaled distance of the hot zone from the top of the potential 
barrier, d = j^. Since the hot zone is in between the top of the barrier and the right min- 
imum, the parameters d and w can only take values between and 1. On the other hand, 
the strength s can take positive values for hot zone and negative values (greater than —1) 
for cold zone. In addition, for each of these cases we have studied the influence of the barrier 
height uq = 7ff- on the equilibration rate. 



The change in the equilibration rate due to the presence of the hot zone is better appreciated 
in terms of a quantity which we call enhancement factor, fb, and defined by 

h = - (18) 
e 

where eo stands for e& when there is no hot zone (s = 0). fb is nothing but the factor by 
which the equilibration rate improves due to the presence of the hot zone. 

An expression for the enhancement factor in terms of all the relevant parameters can be 
found from Eq . ([17|) . By noting that, for a high barrier, k and ki can be expanded about uq 
and Mi, respectively, we find the enhancement factor to be approximately given by 

j. _ exp(sa) cosh(sa) 

1 + s exp(sa) sinh(a) exp(—2u d) 

where a = y^. Even though this equation gives very accurate values comparable to those 
obtained by numerically solving Eq . (fPTD , the above equation can be further approximated 
into a more transparent form when the hot zone is located half way down the top (d = 0.5) 
and is given by 



f _1 



1 + exp ( 2uqw- 



i + (20) 

Equation (|20D shows that for a given barrier height uq and width w, fb increases as the 
strength of the hot zone, s, increases and saturates for large values. In addition, the en- 
hancement factor fb has an exponential dependence on both the width w of the hot zone and 
the barrier height uq. It may be noted that the above expression has the correct limiting 
behaviour as a function of both s — > and w — > 0. Equation (|T9|) shows that the enhance- 
ment factor fb saturates as the blowtorch is moved away from the top of the potential barrier. 

As an independent check of the analytical results (Eq. (|l9|) ), we have numerically solved 
Eq.(|2|) for a given temperature and potential parameters. Before numerically obtaining fb, 
we have studied the evolution of the probability distribution towards the steady state. This 
was done by using an initial probability distribution sharply peaked around the bottom of 
the right well. A 3-dimensional plot of the time evolution of the probability distribution 
towards the steady state is shown in Fig. 2a. The numerically obtained asymptotic steady 
state distribution is shown in Fig. 2b. This clearly has higher population in the left well 
than the right which in the absence of the hot zone would be symmetric about the origin. 
Our numerical results show that this asymptotic distribution is practically identical with the 
steady state distribution P ss (x) given by Eq.@. 

One direct method of determining the long time decay rate of the probability distribution is 
to allow the system to evolve for long enough time towards the steady state. At late stages 
when the probability distribution is observed not to show any appreciable change, we can 
calculate the decay rate by plotting the total probability on the right well as a function of 



time on (In t scale). However, it must be mentioned here that for calculating the equilibra- 
tion rate, a more efficient way is to take the initial distribution to be a distribution which is 
marginally different from the stationary distribution (given by Eq.@). Even so, high values 
of Uq (~ 10) take long computer time for the system to reach the near stationary state re- 
quired for calculating the equilibration rate. For this reason, for the numerical works based 
on Eq.(||) we have limited our calculation to uq = 4. Although this value of uo satisfies the 
high barrier condition, it is on the low side. 

We now consider the influence of the hot zone on the equilibration rate. We shall refer to 
the results obtained from the approximate analytic expression (Eq.([T9|)) as analytical re- 
sults, exact results obtained by numerically solving the root of the equation (Eq. (|17D ) as 
semi-analytical results and the results obtained by numerically solving Eq.(0) as numerical 
results. Figure 3 shows plots of the enhancement factor /& versus strength of the blowtorch s 
for various values of the barrier height uq placing the blowtorch (of width w = 0.1) midway 
from the top of the potential (d = 0.5). (The set of curves a, b and c refer to three different 
values of the barrier height uq = 4, 10 and 15 respectively.) We have shown the analytical 
results by a dashed line, the semi-analytical results by a continuous line and the numerical 
results by filled circles. (For higher values of uq, we have shown only the analytical results 
and semi-analytical results, since numerical results are time consuming.) As can be seen 
from Fig. 3, for higher values of the barrier height Uq, analytical results are not distin- 
guishable from the semi- analytical results over the entire range of s. This clearly reflects 
the fact that large uq approximation made in obtaining Eq.(^) holds good. Consider the 
dependence of /& on the width of the hot zone w. Figure 4 shows plots of /& as a function 
of w for uo = 4, 10 and 15 labelled by (a), (b) and (c) respectively for s — 1 and keeping the 
blowtorch half way down the barrier top, i.e., d = 0.5. As before, for uq = 4, where we have 
numerical results, there is excellent agreement between the numerical results, and analytical 
and semi-analytical results over the entire range of w. For uq = 10 and 15, semi-analytical 
results are indistinguishable from the analytical results. Lastly, we consider the dependence 
of the enhancement factor, /&, on the position of the hot zone from the barrier top, d. Figure 
5 shows plots of fb versus d for uq = 4, 10 and 15 denoted by (a), (b) and (c) respectively 
(keeping the strength of the hot zone s — 1 and its width w = 0.1). For relatively low barrier 
height, Uq = 4, we see that the analytical results agree well with the semi-analytical results 
and the numerical results. Unlike for small barrier height for which /& saturates slowly, for 
large barrier heights ((b) and (c)), /& saturates quickly and stays nearly constant beyond 
d ~ 0.2. Here again, the analytical results agree very well with the semi-analytical results. 
The above results reveal that the equilibration rate is higher when the hot zone is away from 
the barrier top. 

The above dependence of enhancement factor (/&) for the equilibration rate on the position of 
the blowtorch (d), shown in Fig.5, is intriguing. The enhancement factor is seen to saturate 



quickly (to a value greater than unity) as the blowtorch moves away from the barrier top (i.e., 
as d increases beyond while it decreases towards unity as the blowtorch approaches the 
top of the potential barrier. (Of course, for sufficiently small d, the hot zone begins to bracket 
the potential maximum because of finite width.) The above general feature may be physically 



understood in terms of an argument due to Landauer |27j based on the equilibration rate, 
- = — + — , for the occupation of the two wells, where tl and tr refer to the time for 
crossing the barrier from the left and from the right, respectively. Thus, for a blowtorch on 
the right side of the barrier, we expect not only tr to increase as we approach the top of 
the barrier (as d decreases), but also tl should increase as the particles crossing the barrier 
from the left well will be returned back by the blowtorch. This might qualitatively explain 
the result of d dependence of the enhancement factor, However, this argument does not 
enable us to estimate the relative importance of the two effects, namely, the variation of 
Tl and tr with d. To clarify this point, we have carried out calculation of r for a highly 
asymmetric double-well potential in the limit of a very deep right well. We considered two 
cases — (i) when the hot zone is placed to the right of the barrier top and (ii) when it is to 
the left. This enabled us to isolate to a very good approximation the d dependence of tl 
and tr respectively from r(d). Our finding is that the decrease in the enhancement factor, 
fb, as d decreases is dominated by the increase of t r , although r L also increases weakly. This 
is somewhat counterintuitive. 



5 Summary and conclusions 

In summary, we have been able to study the kinetic aspect of Landauer's blowtorch theorem 
using the supersymmetric approach and using a simple model VK-potential. The choice of 
the ^-potential is particularly well adapted to the SUSY method. Our analysis shows that 
the rate of equilibration is substantially improved due to the presence of the hot region. 
The exact magnitude depends on its strength, its width, its location and also on the barrier 
height. We have also obtained an approximate analytical expression for the equilibration 
rate from transfer matrices derived using the SUSY method. These results agree well with 
the results obtained by numerical solution of the associated Smoluchowski equation. 



We expect that the present analysis would be useful in understanding some problems when 
local heating is viewed in a more general context of local noise injection. Viewed from this 
angle, these results are clearly applicable even when fluctuations are athermal. One ex- 
ample where athermal fluctuations play an important role is the depinning of dislocation 
segments from obstacles resulting in movement of dislocations [p8| . Another example where 
the present analysis may be useful is the study of kinetics of phase transformations. In this 
case, the free energy takes the role of the potential and the order parameter takes the role 
of the reaction coordinate. However, it is worth pointing out that in general, subjecting a 
range of order parameter values to excess heating is not easy since there is no correspondence 



between the values taken by the order parameter and its location in space. In this context, 
martensite transformation offers a promising physical situation where we are actually dealing 
with the free energy in terms of strain order parameter. It would be interesting to realize 
the applicability of the present analysis to some experimental situations. 
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Appendix A 



In this appendix we outline the method of evaluating the matrix M that relates the ampli- 
tudes A' 4 , B' 4 of found on the left side of the delta function potential at y = — 1 to the 
amplitudes A , B found on the right side of the delta function potential at y — 1. (See Fig. 
lc.) 



The region of the W potential has four intervals of constant potential V-(x) marked by I, 
II, III and IV as shown in Fig. lc. In these intervals, only simple phase changes in the wave 
function occur. In contrast, the wave function changes in a significant way across regions of 
temperature discontinuity and at points where delta functions are present. First, consider 
the phase changes which are simple to evaluate. The pairs of amplitudes A n+ i and B n+ i 
at the left end of these regions(I to IV) are related to those at the right end, A' n and B' n , 
through the matrix equation 

= P(a)(§) (21) 




• p -ka n 

P(«) = ( n (22) 



where 

T>(n\ = ( e 

e k 

Here a stands for the width of the concerned interval and k is either k or k\ depending on 
whether the interval is outside or within the hot zone. 



Now, consider relating the amplitudes across the points of discontinuity in the temperature 
profile. The probability density P(y, t) which is related to 4>°_{y) through the transformation 
Eq.(Q), is a continuous function of the position including the two points where temperature 
is discontinuous. Therefore, 4>°_{y) on either side of these two points must be discontinuous 
in such a way that P(y,t) remains continuous. In particular, at y\ = d + — , using Eq.(3) 
and using the continuity condition P(y1[,t) = P(yi,t), we get 

<P°Mt) = f^f (23) 

Here y x = y\ ± e with e — > 0. 



The second relation between the amplitudes comes from integrating the SE equation, across 
?/i and exploiting the continuity of P(y,t). This leads to the continuity of {TP). At 
y — 2/1, this implies that 

(0 O (Z/1 + )) , = (0 O -(Z/D) / , (24) 

where the primes denote the spatial derivative. Using Eqs. (|23"D and (0), the pairs of 
amplitudes A%, B\ and A[, B[ found on the two sides of the discontinuous temperature 



profile located at y\ are related through the equation 
where 

J_f(l + s)h + k (l + s^-ko} , , 

1 2k 1 \(l + s)k 1 -k (l + s )k 1 + ko)- 1 ° j 

In the same way, the two sets of amplitudes A2, B2 and A' 2 , B' 2 found on the two sides of 
the discontinuous temperature profile located at y 2 = d — ^ are related through the matrix 
equation 



where 



* 2 r M2 U ) w 



M 1 /fco + (l + s)fci fco-(l + s)A;A , , 

2 2(1 + s)k Uo - (1 + s)h k + (l + s)h ) ■ [Z6) 
We are now left with relating the amplitudes on either sides of the three delta function 
potentials located at y = —1,0 and 1. At all these three delta functions, <f)°_ is continuous 
since temperature is continuous. Thus, the pairs of amplitudes A , B and A' , B' Q on the 
right and the left sides of the delta function located at y — 1 are related through 

A' + B' = A + B . (29) 

Integrating the Euclidean Schrodinger equation across this positive delta function with the 
width of the integration tending to zero gives a second relation between these amplitudes: 

K - B' = (A) + B ) + A - B Q . (30) 

ko 

These two relations gives a matrix equation relating the two sets of amplitudes: 

UH+U)' < 3i » 

The same transfer matrix M + also relates the pairs of amplitudes A' 4 , B' 4 to A 4 , B 4 since 
features around y = 1 are identical to those around y = — 1. When relating the pairs of 
amplitudes A 3 , B 3 to A' 3 , B' 3 corresponding to the two sides of the delta function located 
at y = 0, we note that the only difference is that the sign of the delta function is negative. 
Going through the same procedure as above we get the matrix equation: 



A' 3 
B' 3 




M - n • (32) 



The transfer matrices M± are given by 



M ± 4^ % ,™ ) (33) 



k V ±Uo k ±u . 

Using all the transfer matrices relating the successive sets of amplitudes, the matrix M is 
given by 

in hi 

M = M + P(l)M_P(d )M 2 P(w)M 1 P(l - d )M+. (34) 

2 2 



Figure Captions 



Figure 1: (a) The W - potential, (b) Temperature profile T(x), having the hot zone, (c) 
SUSY partner potential, V-(x). 

Figure 2: (a) Time evolution of the probability distribution P(x, t), towards the steady state 
for a given temperature. The parameter values are : s = 1.0, u$ = 4, w = 0.1 and d = 0.5. 
(b): Very late stage probability distribution obtained by numerically integrating Eq.(||) taken 
as the steady state distribution. Note the dip in the probability value at the location of the 
hot zone shown by an arrow. 

Figure 3: Plot of /& versus s, for d = 0.5, w — 0.1 for the values of (a) u = 4, (b) u = 10 
and (c) m = 15. Dashed lines correspond to the results from Eq . ([19|) , solid lines to the 
semi-analytical method and filled circles to numerical solution. 

Figure 4: Three plots of /& versus w, for d = 0.5, s — 1 for the values of (a) uo = 4, (b) 
uq = 10 and (c) «o = 15. Dashed lines correspond to the results from Eq.([T9|), solid lines to 
the semi-analytical method and filled circles to numerical solution. 

Figure 5: Plots of /& versus d for s = 1.0, w = 0.1 for the values of (a) w = 4, (b) uq = 10 
and (c) u = 15. Dashed lines correspond to the results from Eq.(^), solid lines correspond 
to semi-analytical results and filled circles to numerical solution. 
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